Energy Conversion and Management 78 (2014) 565-573 



ELSEVIER 


Contents lists available at ScienceDirect 

Energy Conversion and Management 

journal homepage: www.elsevier.com/locate/enconman 


Energy 

Conversion 

^Management 




A two control volume model for the Thermal Lag Engine 



CrossMark 


C. Fdez-Aballi Altamirano c ’*, S. Moldenhauer 3 , J. Gonzalez Bayon c , S. Verhelst , M. De Paepe 


a Institute of Thermodynamics and Fluid Mechanics, llmenau University of Technology, P.O. Box 100565, 98684 Ilmenau, Germany 
b Ghent University, Department of Flow, Heat and Combustion Mechanics, Belgium 
c Centro de Estudio de Tecnologias Energeticas Renovables (CETER), Cujae, Cuba 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 6 July 2013 
Accepted 4 November 2013 
Available online 10 December 2013 


Keywords: 

Thermal Lag Engine 
Mathematical model and simulation 


The Thermal Lag Engine (TLE) was proposed by Tailer in 1993 and its mechanical simplicity makes it an 
appealing technology to access renewable sources of energy or recover waste heat. This engine has a sin¬ 
gle moving part, the piston, which delivers power and acts as the prime mover of the working fluid. The 
TLE consists of a hot space with a high heat transfer capacity and a cold space with limited heat transfer 
capabilities. This engine is still not well understood. The paper presents a two control volume model for 
the TLE. In order to validate the model it is parameterized to describe a TLE built by Organ. The simula¬ 
tions of Organ's TLE experiments when compared to his experimental results have a correlation coeffi¬ 
cient of 0.92 and are significantly correlated within the 95% confidence interval. The different 
characteristics of Organ's TLE are discussed based on the simulation results. Conclusions about the 
required characteristics of the heat transfer regimes and energy fluxes for the TLE are made. It is dis¬ 
cussed that the key to a more efficient TLE will result from a careful design of the gas path and the cold 
heat exchanger, so that more of the enthalpy shuttled into the cold space can be converted into work. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Anthropogenic environmental degradation and fuel scarcity 
intensifies the need for economically viable technologies to access 
renewable sources of energy. Since 1807 air engines have provided 
an alternative technology for power generation due to their 
mechanical simplicity and fuel flexibility. James and Robert Stirling 
in 1816 included the regenerator into the air-cycle resulting in the 
patent of the Stirling engine. Other external combustion engines 
following Erickson cycles or Carnot cycles have also been built. 
Yet, it is Stirling technology that has demonstrated high thermody¬ 
namic and mechanical performance proving competitive for niche 
markets such as space applications. The Stirling refrigerator is 
widely used in cryogenic applications [15], 

Stirling cycle prime movers and refrigerators are still an active 
area of research. Researchers are constantly proposing improve¬ 
ments to these machines and better ways to describe them, mostly 
concentrating on the classical configurations: alpha, beta and gam¬ 
ma. The three most challenging practical aspects of Stirling ma¬ 
chines are: (1) lubrication and material selection for the moving 
parts in the hot section, (2) sealing the engine with a moving piston 
and (3) achieving heat exchangers with sufficient heat transfer 
capacity without increasing flow resistance or dead volume. The 
mechanical coupling between the piston and the displacer also 
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introduces mechanical complexities. To address these difficulties 
four different variants of Stirling-like machines have been devel¬ 
oped since the 1960s, the thermoacoustic hybrid Stirling engine 
(TASHE) developed by Swift and his team [12], the free piston Stir¬ 
ling Engine developed by Beale [3], the pulse tube cycle initially 
proposed by Gifford and Longsworth in 1964 [16] and finally the 
Thermal Lag Engine proposed by Tailer in 1993 [13]. These config¬ 
urations reduce the moving parts, eliminating the mechanical cou¬ 
pling between the displacer and the piston, remove the moving 
parts in the hot section of the engine and address the sealing 
problem. 

Tailer was granted a U.S. patent for the TLE in 1995 [14] and 
Fig. 1 illustrates one of his experimental engines. Tailer presents 
an elegantly simple idea for an external combustion engine, which 
uses a single moving part, the piston. The TLE is still not well de¬ 
scribed and understood, yet its mechanical simplicity makes it an 
appealing technology for solar applications or waste heat recovery. 
The purpose of this work is to show that it is possible to describe 
the TLE cycle through a two control volume model of the engine 
and to further the understanding of these machines by studying 
the simulations of the experiments presented by Organ in his book 
[ 11 ]. 

2. Organ and Tailer’s TLE 

In general the TLE consists of two spaces: a hot space where the 
heat to power the engine is supplied and a cold space into which 
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Nomenclature 

P 

density (kg/m 3 ) 

Suffixes refer to 

U 

internal energy (J) 

h hot control volume 

Q 

heat transfer (J) 

c cold control volume 

w 

work (J) 

p piston 

H 

enthalpy (J) 

0 initial or fixed value 

R 

specific gas constant air (J/kg I<) 

max maximum value 

C v 

constant volume heat capacity (J/kg K) 

/ flow path value 

C P 

constant pressure heat capacity (J/kg K) 

me mean value 

h 

heat transfer coefficient (W/m 2 K) 

tot total values of all the volumes 

C c 

inner perimeter of piston chamber 

H and C hot and cold sources characteristics 

Cf 

flow contraction coefficient 

Differential notation f = % 

A p 

Piston area (m 2 ) 

D 

diameter (m) 


L 

length (m) 

Variables and parameters 

d 

thickness (m) 

t time 

X 

Piston displacement (m) 

T temperature (K) 

e 

crank angle starting at top dead center (Rad) 

P pressure (Pa) 

f 

engine frequency (Hz) 

m mass (kg) 

e 

angular velocity (Rad/s) 

V volume (m 3 ) 

r 

Crank radius (m) 

A area (m 2 ) 

l< r 

connecting rod length (m) 


the piston brings the gas intermittently because of its motion. 
Thus, the piston takes the double function of delivering work and 
acting as the prime mover of the fluid in the engine. The TLE piston 
takes the role of the displacer in Stirling devices as its motion in¬ 
duces the movement of the gas through the different sections of 
the engine shifting in time within the cycle when the heat input 
and rejection occurs. Tailer and Organ have both built and mea¬ 
sured TLEs, yet they have different claims about what are the 
mechanisms that generate the net heating and cooling of the gas 
at the right points in the cycle [13,11], 

In Figs. 1 and 2 both Organ and Tailer’s TLE configurations are 
illustrated. In both cases the engine consists of a hot space with 
a large surface area heat exchanger where heat is supplied, fol¬ 
lowed by an interface that separates the hot side form the cold 
side, allowing for the engine material to sustain a temperature dif¬ 
ference. In the case of Organ this is a tube, in the case of Tailer this 
is a flame guard and water jacket. Then the cylinder top is a cooled 
surface, where the hot gas that reaches the piston can dump heat. 
Therefore the arrangement of components in these two engines is 
effectively the same. Yet there are some fundamental differences 
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Fig. 1 . Simplified schematic of Tailer's TLE. 



Fig. 2. Simplified schematic of Organ’s TLE. 


between these machines. The cold heat transfer in Organ’s machine 
is limited by the outside heat transfer coefficient, unlike Tailer’s 
which uses water as the cooling medium. In Organ’s TLE the cold 
heat exchanger is placed before the entry of the gas into the cylin¬ 
der and the bore to stroke length (D p /X max ) ratio is 1.3, Tailer does 
not have a cold heat exchanger there and ( D p /X max ) is 10. Thus dif¬ 
ferent heat rejection behaviors can be expected for these engines. A 
further difference lies in the nature of the interface between heater 
and the cooler, which leads to the two different interpretations of 
how the TLE works. Tailer claims that the TLE should have a small 
as possible interface between the heater and cooler, and that the 
enthalpy shuttling between the two spaces is the important effect 
to consider [13]. Organ states that the TLE should have an optimal 
volume and hydraulic radius interface between the heater and the 
cooler. He refers to this interface as the pulse tube in accordance 
with his claim that the TLE is the inverse cycle of the pulse tube 
cryogenic coolers invented since the 1960s [11], 

Whilst both Tailer and Organ provide a qualitative explanation 
for why their engines work, unfortunately Tailer’s experiments are 
not sufficient to prove his claims and Organ’s simulation does not 
correlate with his experimental results. Therefore the understand¬ 
ing of the thermodynamic phenomena that allows the TLE to pro¬ 
duce work remains uncertain, and will only be unveiled through a 
thorough numerical and experimental investigation as Organ him¬ 
self suggests [11 [.For the Pulse Tube Engine (PTE) such investiga¬ 
tions have been performed by Moldenhauer et al. [5,21], 

3. Attempts to model the TLE 

The previous attempts to model the TLE use a two or greater 
number of control volume description. The earlier attempts by 
Chen and West [8], West [18 and Wicks and Caminero [19] de¬ 
scribe the engine as two volumes, one hot and the other cold. West 
presents the most insightful description of Tailer’s engine and the 
TLE principle in general. In order to do so, he takes simplifying 
assumptions that allow him to identify states and processes within 
the cycle. Thus West presents a quantitative approximation of the 
TLE, indicates the importance of the time dependence of the 
processes in the engine and illustrates the unique effects that come 
from coupling an isothermal gas spring with other processes. 

Wicks and Caminero simulate the TLE by fixing the heat transfer 
regimes inside the engine. They assume nominally constant heat 
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input in the hot section and sinusoidal heat output in the cold sec¬ 
tion. From this they assume that the variation of internal energy in 
the engine also follows a sinusoidal behavior and produce pV and 
TS diagrams. Yet the variation in internal energy is also coupled 
with the work output, and not taking this into account requires 
the use of unrealistic heat transfer regimes to describe the engine, 
as was observed in one of the simulation attempts made, and is 
also referred to in Wicks’s paper [19], 

Organ solves a one-dimensional approximation of the mass, 
momentum and energy balances in the engine. He does not at¬ 
tempt to model Tailer’s engine. Instead he develops a model to 
simulate the engine built and operated by his team, but presents 
no comparison between his simulations and his experimental re¬ 
sults [11], For his simulations Organ assumes laminar flow behav¬ 
iour, and gives special attention to the interface region between 
the cold and hot space of the engine. He refers to this region as 
the pulse tube. He includes a comprehensive description of how 
the pulse tube permits work generation in a constantly heated 
and cooled engine. This is important in Organ’s work because, un¬ 
like Tailer he assumed the cylinder walls and piston head to be adi¬ 
abatic, only transferring heat through the cold stack in the cylinder 
top. Interestingly there is not a pulse tube between the cold and 
hot space in Tailer’s engine. Nevertheless Organ claims that the 
TLE corresponds to the pulse tube engines and refrigerators family 
already patented since 1964 [11], 

The concept of the Pulse Tube Engine (PTE) presented by Ham- 
aguchi et al. [4] has been studied both numerically and experimen¬ 
tally. Yoshida et al. [20] in 2009 measured the work flux in a pulse 
tube engine showing they could be classified into the thermoacou¬ 
stic standing wave engine group. In 2012 Moldenhauer et al. pre¬ 
sented thorough numerical and experimental studies of the PTE 
reporting measured indicated powers and efficiencies of the order 
of 6W and 8% [6,5], These results reaffirm the PTE’s potential for 
low grade heat recovery applications. 

The work presented here returns to West’s simulation approach 
by simplifying the machine into two control volumes. As in West’s 
and Wicks’s work one is a hot volume with high heat transfer 
capacity and the other is a cold volume with limited heat transfer 
capacity and varying heat transfer area. Our modelling approach, 
unlike the previous attempts does not pre-establish the type of 
thermodynamic processes making up the cycle. On the other hand 
a dynamical modelling approach is taken; integrating energy, 
momentum and mass balances between the volumes in time. This 
model is the basis of the simulations later presented. 

4. The model 

In this model the engine volume consists of two spaces, one 
with fixed dimensions exposed to the hot source and one with 
varying dimensions exposed to the cold source. The total dimen¬ 
sions of the cold space’s volume and heat exchange surface are 
determined by the piston’s displacement. The displacement of 
the piston is as if it was coupled to a crankshaft and flywheel 
mechanism with constant angular velocity. 

It is the general consensus that unlike thermoacoustic ma¬ 
chines, the significant effects in the TLE are derived from bulk fluid 
motion [2], Then it should be possible to describe the engine with 
the interaction of two control volumes with different heat transfer 
regimes coupled by an enthalpy flux. This is the general treatment 
given to the TLE models found in literature, with the exception of 
Organ who attempts to capture the effects arising from the tem¬ 
perature profiles in the gas and engine walls [11], 

The temperatures of the hot source and cold sink are constant. 
The working fluid is air and is assumed to behave like an ideal gas. 
In most of the TLE engines measured and modelled that have been 
found in literature, the temperatures of operation of the gas are 


such that the ideal gas behaviour is an adequate approximation. 
For air at the range of operation temperature and pressure of the 
engines there is less than 10% variation in c p and c„ [9], thus allow¬ 
ing for these values to be treated as constants. The same approxi¬ 
mations are taken by other authors. 

The schematic of the model is presented in Fig. 3. 

The TLE can be simplified in three main spaces, a heated space a 
with a very large heat transfer capacity, a cooled space with vary¬ 
ing volume and limited heat transfer capacity and an adiabatic 
space connecting the heater with the cooled space. The gas in the 
heater, because of its heat transfer capabilities, remains always 
close to the source temperature, whilst the gas leaving the heater 
into the adiabatic space takes a different thermal behaviour deter¬ 
mined by the enthalpy flux from the heater, the change in volume 
of the engine and the heat rejection in the cooled space. Thus the 
volume of the adiabatic interface can be added into the cold 
volume. 

In the model the volume in the cold CV is varied by the motion 
of the piston, which also varies the cold heat transfer area. The two 
volumes are connected by a flow area, through which enthalpy is 
shuttled between them. Mass, energy and momentum balances 
are applied to the control volumes performing a dynamical simula¬ 
tion. The internal energy of the gas in the volumes cycles in time, 
as the gas is brought out of equilibrium with its surroundings by 
the piston’s motion. Fig. 3 illustrates the model’s arrangement 
and from it the energy balances of the hot and cold volumes can 
be readily stated. The convention is that the energy fluxes leaving 
the volumes are taken to be negative. Then the energy balances for 
the CVs are: 

For the hot volume: 

dU/dt = Qh + Hh (1) 

For the cold volume: 

dU c /dt = Q c - Hh - W (2) 

The dependence on the temperature gradient and the heat 
transfer coefficients are both of importance for the model to esti¬ 
mate the heat input and rejection. The exterior of the engine and 
the wall materials can be readily selected and modified to increase 
heat transfer capacity. Air cooled engines usually have their cool¬ 
ing capacity limited by the heat transfer coefficient of the exterior 
air. In the cases where the heat transfer outside is enhanced, the 
limitation to heat transfer will come from the internal resistances. 
Other authors have expressed the importance of estimating suit¬ 
ably the values of the heat transfer [11,19], 

The heat exchange in the hot volume can be expressed as: 

(b=A h h h (T H -T h ) (3) 

For the cold volume there is a fixed heat exchange area and a 
variable one, which is varied by the Piston’s motion. The cold fixed 
area is usually composed of the cylinder head, the piston surface 
and in some cases a heat exchanger included after the interface 
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region and before the piston chamber. Examples of inclusions that 
increase the fixed cold heat transfer area have been mentioned 
both by Tailer and Organ [14,11], Even though within one control 
volume the heat transfer phenomena is averaged in space, special 
attention needs to be placed on what heat transfer mechanism is 
taking place in what surfaces, such that the magnitude of the heat 
transfer occurring can be approximated correctly. The piston, even 
if it is receiving an impinging velocity, is usually made of a poor 
conducting material and the heat transfer is limited by the outside 
heat transfer coefficient. Thus the piston area can be neglected or 
taken to have a lower heat transfer capacity. Further, in surfaces 
perpendicular to the predominant flow direction, because of flow 
separation generated by entry and exit effects the heat transfer is 
also limited. Therefore in the model some of these areas can be ne¬ 
glected or should have their heat transfer capacity limited with re¬ 
spect to other areas. These effects are introduced through a 
correction factor a 0 taking values between 0 and 1. The value of 
a 0 should be 0 when the heat transfer can be neglected, 1 to con¬ 
sider the complete heat transfer capacity and a middle value when 
the combined effect of the different heat transfer characteristics of 
these areas needs to be considered. 

The inclusions before the piston tend to be flow through heat 
exchangers with a heat transfer area A xc , and have a significant 
contribution to the heat rejection in the cycle. The usual position 
of this heat exchanger after the adiabatic space and before the 
cooled piston chamber makes it important to consider the temper¬ 
ature of the gas entering it. During the expansion this heat exchan¬ 
ger is receiving hot gas coming from the heater after flowing 
through the adiabatic space and in the compression it is receiving 
cold gas from the piston chamber. In order to describe this effect an 
extra term can be included into the cold heat transfer expression, 
where the temperature of the gas entering the heat exchanger cor¬ 
responds to the hot temperature or when the adiabatic space exists 
to the mean temperature value for the whole engine T, 0 , during the 
expansion and to the cold volume temperature value T c during the 
compression. 

The piston motion exposes the chamber, that has a circumfer¬ 
ence C c , such that the piston’s displacement multiplied by the cir¬ 
cumference gives the variable heat exchange area. In piston 
configurations where the bore D p to stroke length X max ratio is close 
to 1 or smaller, this variable lateral area contributes to the heat 
transfer dynamics that allows for work production in the engine, 
creating maximum heat rejection capacity at maximum engine 
volume. In the cases where D p /X max is much greater than 1 the flow 
in the piston takes a large radial component and other phenomena 
become the important effect, the heat transfer model presented 
here does not consider these cases. Then the heat transfer in the 
cold volume can be expressed by: 

a f (Moc + C c x)h c (T c - T c ) +A xc h c (T c - T tot ) if x> 0 

^ l (Moc + C c x)h c (Tc - T c ) + A xc h c (T c -Tc) if x< 0 

The value of T cot is calculated by: 


Uh + Uc 

CvOAtot 


( 5 ) 


The cold and hot internal heat transferred coefficients can be 
modelled using correlations found in literature taking time aver¬ 
aged mean values for the gas properties between the hot and cold 
CV. For the TLE the hot CV heat exchanger usually is a porous med¬ 
ium. There is no flow through the volume and thus the contribu¬ 
tion of bulk flow velocities to the heat transfer phenomena can 
be neglected. Then the predominant heat transfer mechanism 
can be assumed to be natural convection in a medium of a small 
hydraulic radius. To estimate the coefficient value the Churchill 
and Chu correlation is used [17], where the characteristic length 


for the correlation corresponds to the hydraulic diameter D h of 
the porous medium of the hot heat exchanger. For the cold heat ex¬ 
change coefficient the Annand correlation is used [1 . This correla¬ 
tion is specific to heat transfer mechanisms in cylinder-piston 
arrangements. In the TLE only gas compression and expansion pro¬ 
cesses take place, then the radiation term in the correlation can be 
neglected, and for initial approximations the correction factor a in 
the correlation can be taken to be the mean value of 0.5 [ 1 ]. 

Taking air standard analysis, the enthalpy fluxes can be 
expressed by: 

H h = c„m h T h/c (6) 

Because mass is conserved between both volumes, the enthalpy 
flux term in the cold energy balance is just the opposite of the hot 
energy balance. The T h/C variable in Eq. (6) corresponds to the 
temperature of the CV from which the mass flux is coming. In 
the cases when the effect of the interface region needs to be con¬ 
sidered and a better approximation is achieved when T h/C takes 
the value of T tot instead of T c during the compression. The direction 
of the mass flux is determined by the pressure difference between 
the volumes. 

A mass flux model can be obtained by treating the flow between 
the volumes as if it was undergoing a sudden contraction or expan¬ 
sion between two spaces at different pressures. The mass flux 
model used to describe the mass transfer between the volumes 
in the TLE is derived from the integrated form of the momentum 
equation, which by making the following assumptions yields Eq. 
(7), where A 2 represents an area downstream of the flow from Ai 
[17]: 

1. The flow is considered to have a mean density for each time 
step. 

2. The flow is steady between time steps. 

3. The gas is always parting from an initial velocity of zero. 

4. The flow is one dimensional. 

5. The mass flux between the volumes will be modelled over 
the interconnecting flow area. 

6. No disipative losses are considered. They will be incorpo¬ 
rated later using a loss coefficient C f in the mass flux 
equation. 

7. The flow has no change in gravitational potential energy. 

A 2 {P,-P 2 ) = p me {A 2 v 2 2 -A,v\) (7) 

From Eq. (7) a model can be derived to calculate the flow velocity 
for every time step proportional to the pressure differences be¬ 
tween the hot and cold CVs. The value of the flow density corre¬ 
sponds to a linearisation of the density in space, performed by 
taking the mean density value between the volumes. Finally as 
pressure differences between the volumes range around 100 Pa, 
the flow is reciprocating and there is no flow through any of the vol¬ 
umes, it is possible to assume the initial velocity of every flow pro¬ 
cess to be zero. Changing to the appropriate nomenclature and 
taking the zero initial velocity assumption, Eq. (7) can be simplified 
to get the flow velocity between the hot and cold volumes. In the 
reciprocating case where the flow is changing direction as the pres¬ 
sure difference between the control volumes inverts, A 2 becomes A f 
given that in one direction or the other, the flow is parting from the 
volume with the highest pressure, towards the flow area separating 
the two spaces. Ai becomes A 0 because it would be the cross- 
sectional area of the hot or cold CV depending of the direction of 
the pressure gradient. As the initial flow velocity is assumed to be 
zero, A 0 disappears from the equation making the flow velocity only 
dependent on the flow area between the control volumes. 

A f (P h -P c ) = p me (A f vj-A 0 v 2 0 ) (8) 
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vl = o 


Vf- 


(Ph-Pc) 


To get to the mass flux equation this velocity needs to be mul¬ 
tiplied by the flow area A f for obtaining a volumetric flow rate and 
then by the mean density p me . As the flow between the volumes 
can be occurring through some interface passage of smaller diam¬ 
eter than the piston bore, the gas is flowing into the piston cham¬ 
ber through a sudden expansion and out of it through a sudden 
contraction. Empirical correlations for sudden expansion and con¬ 
tractions introduce a correction factor C f corresponding to the 
losses originated by the entry geometry. For the sudden expansion 
this value is taken to be 1 and for the sudden contraction 0.6 [7], by 
moving the value of this coefficient closer to zero further dissipa¬ 
tive losses can be introduced into the model. Further, it has been 
shown experimentally that to use a mean density value between 
the two spaces for the flow density is an adequate approximation 
[7], This then results in the expression for mass transfer used in 
the model, where the flux leaving the volume is considered 
negative. 


m h 


(-Cf hc p me A f] /&^ifP h -P c >0 

\c fchPme A fS /^ifP h -Pc<0 


( 9 ) 


This discrete equation is required because both the flow direc¬ 
tion and the boundary conditions vary discretely with the direction 
of the pressure gradient. The values of P c and P h are also swapped 
between one instance of the equation and the other, thus the solu¬ 
tion of the square root is always defined. 

To complete the enthalpy transfer expression a discrete 
relationship is needed defining the temperature of the mass 
that is transferred and thus its enthalpy. This is expressed in 
Eq. (10). 


From the definition of work, the shaft power term in Eq. (13) can be 
expressed in terms of the pressure in the cold chamber, the piston’s 
area and the piston’s velocity. 

I (a 0 Ao c + C c x)h c (Tc-T c )+A xc h c (Tc-T tot )-c p mi,T h/c -P c A p x ifx^O 
\ (,a 0 Aoc + C c x)h c (T c - T c )+A xc h c (T c - T c ) - c p m h T h/c - P c A p x if x < 0 

(14) 


dU c /dt 


The model for the piston motion corresponds to the standard 
equations of a crankshaft and flywheel mechanism where 
the flywheel remains at a constant angular velocity: 


x=k + r- \/ k 2 - r 2 sin 2 8 

Differentiating Eq. (15) results in: 


x = rsindd + 


8r 2 sin8cosd 
\Jk 2 - r 2 sin 2 8 


(15) 


(16) 


corresponds to the given engine fre- 


The angular velocity 
quency/in Hz: 

8=f2n (17) 

The system consists of 10 variables. It can be solved by adding 
the perfect gas law, the definition of the internal energy and an 
expression for the conservation of mass, to Eqs. 9, 10, 12, 14, and 
16. Through simple algebraic manipulations the pressure and the 
temperature can be expressed by the internal energy of the vol¬ 
umes. This procedure allows simplifying the integration to three 
variables: U h ,U c ,m h . This guarantees mass conservation and re¬ 
duces numerical error in the integration scheme. Eqs. (12) and 
(14) become: 


dU h /dt = A h h h (T H - -^-) + c„m h T h/ , 
c v m h 


(18) 


dU /dt={ 

\ (aoAoc + C c x)h c (T c -^)+A xc h c (T c ~^)-c P m h T h/c -^y^A p xif x<0 


dTh/c/dt 


T h if P h - P c > 0 
T c if P h - P c < 0 


( 10 ) 


In the cases where the interface exists Eq. (10) becomes: 



T h if P h -P c > 0 

T t ot if Ph — Pc < 0 


( 11 ) 


Then Eqs. 9, 18 and 19 compose the differential system to be 
integrated. As the angular velocity of the flywheel is assumed con¬ 
stant, the velocity of the piston and the cold volume could be read¬ 
ily solved from Eq. (15)-(17). 

5. The integration scheme 


These conditional statements add numerical instabilities at 
the point when the pressure gradient shifts direction. This hap¬ 
pens because at this point an abrupt change occurs in the first 
derivative of the mass flux, which is the second order derivative 
of the CV masses. On the other hand, they describe effectively 
the fact that hot mass comes into the cold volume from the 
heater and cold mass enters the heater from the interface or 
the piston chamber directly. For example in Tailer’s engine, Eq. 
(10) is adequate, as the interface region between the hot and 
cold CV is not appreciable. Considering Organ’s engine with a 
large interface between the heater and the cooler Eq. (11) is 
more appropriate. 

Then, substituting Eqs. (3) and (6) in Eq. (1) as well as Eqs. (4) 
and (6) in Eq. (2) we have the following expressions for the internal 
energy of the hot and cold volume. 

dU h /dt = A h h h (T H -Tf) + c p m h T h/c (12) 

(a 0 A 0c + C c x)h c {T c - T c ) +A xc h c (T c - T tot ) - c p mj h/c -Wifx^O 
(Moc + C c x)h c {T c - T c ) +A xc h c (T c - T c ) - c p m h T h/c - W if x < 0 

(13) 


The system is integrated with a variable order solver specifically 
used for solving stiff differential equations. The integration vari¬ 
able is time, and the simulation code is designed to run until a stea¬ 
dy state is reached, having a residual error between cycles below 
1% of the indicated power calculated. The energy balances of the 
cycles are also checked to have a divergence smaller than this va¬ 
lue. The code has been designed to be able to iterate through differ¬ 
ent parameter values and calculate the cycle indicators making it 
possible to study the influence of different parameters in the func¬ 
tioning of the engine. 

6. Modeling Organ’s experiments 

The parametrization of Organ’s TLE for the two control volume 
model is presented next. A comparison of the simulation results to 
Organ’s experiment is presented in Fig. 5. The comparison demon¬ 
strates the validity of the modelling approach. 

Fig. 4 shows the schematic of Organ’s TLE and what regions of 
the engine are considered to be inside the hot and the cold control 
volumes. To simulate using two control volumes the pulse tube 
space was added to the cold control volume and the pulse tube 
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Fig. 4. Control volumes for Organ’s TLE. 


walls were considered to be adiabatic. The cold heat exchanger 
space was also added as a fixed volume to the cold CV. The varying 
cold heat transfer area in the piston chamber was also considered 
in the simulation in contrast with Organ’s assumption of an adia¬ 
batic cylinder [11], Due to the poor conductivity of the piston 
material the heat transfer through that surface was neglected. Con¬ 
sidering the size of the interface between the heater and piston Eq. 
(11) was used to describe the enthalpy flux between the volumes. 
The parametrization used for Organ’s engine in the simulation can 
be seen in Table 1. 

In the case of Organ’s TLE the largest resistance to the heat 
rejection is on the ambient side of the cold heat exchange surfaces. 
The external heat transfer coefficient was estimated to be 20 
W/m 2 K, a standard value for a poor force convection or good free 


convection coefficient for air. To introduce the effect of the external 
heat transfer coefficient the cold heat transfer model was changed 
to Eq. (20). 

i 1 , i' ( a oA)c + C c x)(Tc — T c ) + ■■—i. 1 —-i- (T( — T tot ) if x ^ 0 

kout TC >WW h c A xc 

(a 0 A 0c + C c x)(T c - T c ) + , 1 , , (T c - T c ) if x < 0 

kout Tic h cout A cout h c A xc 

( 20 ) 

In his book Organ presents measurements of the shaft power of his 
TLE keeping a constant electrical power input and changing the en¬ 
gine’s operation frequency. In order to describe the experiment a 
model for mechanical losses was introduced in the simulation using 
the standard analysis for slider crank mechanisms [10]. A compari¬ 
son of the results of the model with Organ’s experiments can be 
seen in Fig. 5. The experimental data has a correlation coefficient 
of 0.92 with the model, yielding a statistically significant correlation 
within the 95% confidence interval. This shows that the proposed 
modelling approach can be adequate to describe Organ’s TLE 
experiments. 


Table 1 

Parameterization of Organ’s TLE. 


Parameter 


Dimensions 

Value 

Hot control volume 

V „ 

Hot volume 

m 3 

9.30 x 10~ 6 

Pe 

Electrical input power 

W 

4.30 x 10 1 

Th 

Hot source temperature 

K 

9.00 x 10 2 

Ah 

Hot heat exchange area 

m 2 

4.49 x 10- 1 

hh 

Hot heat transfer coefficient at 9.5 Hz 

W/m 2 1< 

3.14 x 10 2 

D„ 

Hot heat exchanger hydraulic diameter 

m 

8.00 x 10~ 5 

k, 

Conductivity of the insulation material 

W/m 

4.00 x 10~ 2 

Th, 

Thickness of the insulation material 

m 

2.50 x 10~ 2 

Cold control volume 

Xo 

Minimum displacement 

m 

1.00 x 10~ 3 

Xmax 

stroke length 

m 

2.47 x 10~ 2 

A p 

Piston head area 

m 2 

8.04 x 10~ 4 

Voc 

Minimum cold volume 

m 3 

1.14 x 10~ 5 

Tc 

Cold sink temperature 

K 

3.00 x 10 2 

Aq c 

Piston chamber fixed heat transfer area 

m 2 

8.55 x 10~ 4 

A xc 

Cold heat exchanger area 

m 2 

1.62 x 10~ 3 

Axcout 

Exterior cold heat exchanger area 

m 2 

1.55 x 10~ 2 

Cc 

Cylinder circumference 

m 

1.02 x 10~' 

he 

Cold heat transfer coefficient at 9.5 Hz 

W/m 2 K 

7.02 x 10 1 

hcout 

Ambient heat transfer coefficient 

W/m 2 1< 

2.00 x 10 1 

a 

Annand correlation coefficient 


5.00 x 10“’ 

flo 

Fixed area correction factor 


1.00 x 10° 

Flow passage between the 

Cfhc 

volumes 

Flow coefficient from hot to cold 


1.00 x 10° 

Cfch 

Flow coefficient from cold to hot 


1.00 x 10° 

A f 

Flow area 

m 2 

5.03 x 10~ 5 

Mechanical parameters 

f 

Engine frequency 

Hz 

9.50 x 10° 

k 

Connecting rod length 

m 

6.00 x 10~ 2 

r 

Crank length 

m 

1.24 x 10~ 2 

c tp 

Piston friction coefficient 


2.00 x 10~' 

C J5 

Bearings friction coefficient 


5.00 x 10~ 4 

General parameters 

Tfni 

Filling temperature at BDC 

I< 

3.00 x 10 2 

Pm 

Filling pressure at BDC 

Pa 

1.00 x 10 5 

Cp 

Specific heat capacity at constant pressure 

J/kg K 

1.01 x 10 3 

Cy 

Specific heat capacity at constant volume 

J/kg K 

7.18 x 10 2 

R 

Gas constant 

J/kg I< 

2.87 x 10 2 
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Fig. 5. Comparison of the measurements of Organ’s Thermal Lag Engine and the 
two control volume model simulation [11]. 

7. Investigating Organ’s TLE cycle 

The model presented above can also be used to investigate the 
cycle behaviour of Organ’s TLE. Cycle simulations have been per¬ 
formed using the previous model and parametrization. The engine 
was simulated at frequency values of 7 Hz, 9.5 Hz and 17.5 Hz cor¬ 
responding to points in the experiments. The cycle indicators were 
calculated and the cyclic behaviour of the thermodynamic proper¬ 
ties studied. From the description achieved with a two control vol¬ 
ume model it is explained what factors can be considered for 
improving the performance of TLEs. Table 2 shows the cycle indi¬ 
cators for the different frequencies. 

From Fig. 5 and Table 2 it can be seen that the engine perfor¬ 
mance is 5 times less with a 2-fold increase in frequency. The inter¬ 
nal heat transfer resistance is reduced with increased piston 
velocity, yet the outside resistance remains unaltered limiting 
the heat rejection capability. Therefore, while the cycle period is 
reduced by a factor of 3 the heat transfer capacity is intensified 
only by 50%, resulting in an increase of the mean cycle temperature 
with engine frequency, as it can be seen in Fig. 7. Then the inten¬ 
sification in heat transfer observed at higher frequencies is mainly 
due to the increase in mean engine temperature. In this engine the 
limited heat rejection capacity is one of the causes for the poor per¬ 
formance. This can be appreciated in the indicated diagrams in 
Fig. 6. The hot CV has sufficient heat transfer capacity maintaining 
the same gradient for the expansion process at all frequencies. Yet 
the compression curves cross the expansion lines earlier in the cy¬ 
cle reducing work output as the engine speeds up. This effect is 


x 10 5 



Engine volume (cc) 

Fig. 6. Pressure vs. volume diagrams for Organ's TLE at different engine frequencies. 


generated by the reduction in cycle period and the limited cold 
heat transfer capacity as West indicated in his theoretical work 

[18], 

Therefore the external limit to heat rejection hinders engine 
performance. The effect of insufficient cooling can be best observed 
in Figs. 8 and 9, where the heat transfer and temperature dynamics 
in the hot control volume is presented. There it can be seen that 
heat delivery capacity in this volume allows the gas to maintain 
near isothermal behaviour as predicted by West and Chen 8,18], 
Fig. 9 shows that when the piston reaches bottom dead center 
(BDC) the gas in the heater has reached the heater temperature. 
Then when heat rejection is insufficient, and as the heater delivers 
heat to the incoming cold gas during the compression, the temper¬ 
ature in the hot CV raises above the source temperature. For this 
part of the cycle the engine is pumping heat, attempting against 
power delivery. Using a regeneration strategy to recover this en¬ 
ergy is inadequate as it does not make thermodynamic sense to 
convert work into heat for later recovery. Then as the heater with 
a large heat transfer capacity is a desired attribute in the engine, it 
seems a much better strategy to design the heat rejection mecha¬ 
nism such that this effect does not take place. Thus generating the 
triangular pV diagrams that West [18] predicted would maximize 
power delivery in a TLE. 


Table 2 

Cycle and engine indicators for Organ's TLE. 


Indicators 

Units 

Values 



Engine frequency (f ) 

Hz 

7 

9.5 

17.5 

Heat transfer in hot CV (Q h ) 

j 

2.95 

2.81 

2.43 

Heat transfer in cold CV (Q c ) 

j 

-2.89 

-2.76 

-2.42 

AH from the hot volume 

j 

-2.95 

-2.80 

-2.43 

Indicated work (W,) 

j 

6.39 x icr 2 

5.37 X 10~ 2 

9.60 x 10- 3 

Hot heat transfer coefficient ( h h ) 

W/m 2 K 

314 

314 

314 

Cold heat transfer coefficient (h c ) 

W/m 2 K 

56.7 

70.3 

108.0 

Hot CV heat transfer capacity 

W/K 

141 

141 

141 

Cold heat exchanger heat transfer capacity 

W/K 

7.09 x 10~ 2 

8.32 x 10~ 2 

1.12 X 10-’ 

Piston chamber mean heat transfer capacity 

W/I< 

3.21 x 10~ 2 

3.32 x 10~ 2 

3.64 x 10~ 2 

Indicated power (P,) 

w 

4.48 x 10~' 

5.10 x 10-' 

1.68 x 10-' 

Shaft power (P s ) 

w 

3.88 x 10~' 

4.29 x 10-' 

3.29 x 10~ 2 

Indicated efficiency (77 = Wj/Qj,) 


2.17% 

1.91% 

0.40% 

Volume per watt ( Vtot/Pi ) 

cm 3 /W 

91.5 

79.4 

241.0 

A utot error 

J 

1.50 x 10~ 3 

1.90 x 10- 3 

2.10 x 10“ 3 
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Fig. 7. Total temperature for Organ's TLE at different engine frequencies. 



Crank angle (Rad) 

Fig. 8. Heat transfer in the hot control volume for Organ’s TLE at different engine 
frequencies. 



Fig. 9. Temperature in the hot control volume for Organ’s TLE at different engine 
frequencies. 


Simply dumping sufficient heat will improve engine perfor¬ 
mance. In general the limit to heat transfer in the engine is created 
by the external coefficient, yet enhancing external heat transfer 
can be readily done. Internally enhancing the heat transfer in the 
engine, especially in the piston chamber is difficult, postulating 
the biggest limitation for sufficient heat rejection. As in Organ, 
many TLE’s include a cold heat exchanger before the piston where 
it is easier to enhance heat rejection. Yet in engines it is not only 
important how much heat is being dumped, but also when in the 
cycle the heat rejection occurs. When a heat exchanger is included 
before the piston, heat is rejected before it can be transformed into 
work hindering power delivery. For best performance maximum 
heat rejection should occur towards BDC. Even though Organs 
TLE is net heating towards top dead center (TDC) and net cooling 
towards bottom dead center (BDC), the cooling capacity in the 
gas path before the cylinder, governed by the cycle temperature 
dynamics, shifts maximum heat rejection towards TDC as seen in 
Fig. 10 and Fig. 11. 

The cold heat transfer dynamics observed in Fig. 10 has the 
undesired effect of dumping the enthalpy delivered from the hea¬ 
ter before it can be converted into work. This effect happens more 
intensely as the engine speeds up given that hotter gas is 



Crank angle (Rad) 


Fig. 10. Heat transfer in the cold control volume for Organ's TLE at different engine 
frequencies. 



Fig. 11. Temperature in the cold control volume for Organ’s TLE with different 
engine frequencies. 
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transported into the cold heat exchanger. This can be improved by 
eliminating the cold heat exchanger from the gas path before the 
heater. Then sufficient heat rejection can be attempted by enhanc¬ 
ing heat transfer in the piston chamber. Considering that the vari¬ 
able exposure of the cold area in that space will shift maximum 
cooling towards BDC this is an ideal improvement. This is the ori¬ 
ginal idea suggested by Tailer and developed by West [13,18], 

Also the gas temperature dynamics in the hot CV is dominated 
by the heat transfer capacity of the heater, yet outside this space, 
the gas temperature dynamics is determined by the volume 
change. The cold heat exchangers always receives gas at lower 
temperatures towards BDC further reducing the heat transfer capa¬ 
bilities at that moment. If the cold heat exchanger was placed di¬ 
rectly after the heater to change this, it would receive the hottest 
gas at the beginning of the expansion, thus simply incrementing 
the magnitude of the cooling at TDC and reducing even further 
the work generation. A way to improve that has been to separate 
the heater from the cold heat exchanger, for instance using the 
interface space. It is possible to optimize this space to achieve 
the right cyclic magnitude of heat rejection thus already improving 
engine performance, yet it might be difficult to use this strategy to 
shift in time when the maximum heat rejection occurs. Thus, in 
consequence with the removal of the cold heat exchanger from 
the gas entry into the piston chamber, it might be reasonable to 
eliminate the interface all together. This would allow reducing 
dead volume in the engine. Optimization of the heat transfer in 
the piston chamber could yield more satisfactory results. 

Enhancing the heat transfer in the piston chamber to take 
advantage of the area variation is difficult. One idea could be to in¬ 
clude micro-fins towards the BDC section of the chamber, guaran¬ 
teeing sufficient heat transfer at right time in the cycle. Other 
authors have suggested different ways of achieving maximum heat 
rejection closer to BDC, including open ports [8 and optimizing the 
interface between the heater and cooler [11]. In the pulse tube en¬ 
gine variant other ideas to improve performance have been, to 
optimize the compression ratio, change the filling pressure and 
to include a bypass between the heater and the cold heat exchan¬ 
ger before the piston [5,4,20], In the TLE a bypass could provide a 
parallel enthalpy path that does not pass through the cold heat ex¬ 
changer delivering more energy directly for the expansion. Using a 
strategy like this may allow for optimizing the flow resistances 
such that only the desired amount of energy is dumped and more 
of the enthalpy shuttled between the volumes gets converted into 
work. 

8. Conclusions 

A two control volume model for the Thermal Lag Engine has 
been postulated. The work demonstrates that it is possible to han¬ 
dle Organ’s pulse tube interface and simulate his experimental 
results with this approach. The simulations of Organ’s TLE experi¬ 
ments, when compared to with his experimental results, have a 
correlation coefficient of 0.92 and are significantly correlated with¬ 
in the 95% confidence interval. 

This model is used to simulate and study the cycle behaviour of 
Organ’s TLE. It is concluded that while Organ’s TLE has an adequate 
heating capacity, it has insufficient cooling capabilities limiting 
work production at higher frequencies. Also it is discussed that 
the heat rejection in this TLE happens at TDC and not when it is 
best at BDC. 

The discussion shows that for the engine to perform well the 
hot heat exchanger should have the largest possible heat transfer 
capabilities, whilst the cold heat transfer should be limited to an 
optimal value and in time within the cycle. 


Also it is concluded that the regeneration in the hot control vol¬ 
ume, due to insufficient cooling, is an undesired effect and should 
be avoided. 

The analysis gives that it would be beneficial, in order to in¬ 
crease power output and efficiency, to consider the gas path within 
the engine in order to shift the cooling capacity toward BDC. It is 
shown how the cold heat exchanger placed in the gas path before 
the piston does not attain this effect and hinders the performance 
of the engine, even if it allows guaranteeing the required magni¬ 
tude of the heat rejection. Thus it is discussed that it might be ben¬ 
eficial to eliminate the cold heat exchanger before the piston 
completely and with it the pulse tube interface also reducing dead 
space. It is possible to shift the cold heat transfer to happen at BDC 
by having the piston’s motion dominate the magnitude of the 
available heat transfer area. This allows concluding that the engine 
performance can be improved by enhancing and optimizing the 
heat transfer in the piston chamber. 
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